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Abstract. Some questions raised by _fermj-LAT data about blazars are summarized, along 
with attempts at solutions within the context of leptonic models. These include both spectral 
and statistical questions, including the origin of the GeV breaks in low-synchrotron peaked 
blazars, the location of the gamma-ray emission sites, the correlations in the spectral energy 
distributions with luminosity, and the difficulty of synchrotron/SSC models to fit the spectra 
of some TeV blazars. 



1. Introduction 

The 11 June 2008 launch of the Gamma-ray Large Area Space Telescope, subsequently renamed 
the Fermi Gamma ray Space Telescope, opened a new era into studies of radio galaxies and 
blazars. The Large Area Telescope (LAT) on Fermi, with a peak on-axis effective of ~ 8000 cm^ 
and a single photon angular resolution for 68% containment at 1 GeV of ~ 0.8° [1], represents 
an improvement by 1 - 2 orders of magnitude in sensitivity over its predecessors EGRET on the 
Compton Gamma-ray Observatory and AGILE [2]. In the Third EGRET Catalog of Gamma 
Ray Sources [3j , 66 high-confidence ( ^ 5a) detections of blazars and one probable radio galaxy, 
Centaurus A were reported. Ten high-latitude sources at high galactic latitude (|6| > 10°) were 
detected with significance > lOcr. Now, just over three years into the Fermi mission, three 
important lists of Fermi-hAT AGNs have been or are about to be published: 

(i) The LAT Bright AGN Sample (LBAS) [3]. Based on the first three months of scientific 
operations of the Fermi satellite, the LBAS consists of 106 high-latitude sources with 
TS ^ 100, or significance ^ lOcj. It consists of 58 flat spectrum radio quasars (FSRQs), 
42 BL Lac objects, 2 radio galaxies, namely, Centaurus A and NGC 1275 (3C 84), and 4 
blazars with unknown classification. 

(ii) The First LAT AGN catalog (ILAC) [5]. This catalog contains 671 sources that 
are associated at high confidence with 709 AGNs (some 7-ray sources have multiple 
associations). The associations are based on counterpart catalogs, such as the flat-spectrum 
8.4 GHz CRATES (Combined Radio All-Sky Targeted Eight GHz Survey) Catalog (6] and 
the Roma BZCAT blazar catalog [7J. The ILAC clean sample is a subset of the high- 
confidence associations, and consists of 599 AGNs with no multiple associations or other 
analysis flags. These sources subdivide into 275 BL Lac objects, 248 FSRQs, 50 blazars of 
unknown type (either because a source lacks an optical spectrum, or the optical spectrum 



is inadequate to determine if it is a BL Lac or FSRQ), and 26 other AGNs. These "other" 
AGNs include 7 radio galaxies/misaligned AGNs, 4 radio-loud narrow line Seyfert 1 (RL- 
NLSyl) galaxies, 4 narrow line radio galaxies, 3 starburst galaxies (which of course are not 
actually AGNs), and 9 other sources. 

(iii) The Second LAT AGN catalog (2LAC) [8j. The 2LAC contains 1016 sources that are 
associated at high confidence with AGNs. The 2LAC clean sample consists of 885 sources 
(46% larger than the ILAC clean sample), consisting of 395 BL Lacs, 310 FSRQs, 156 
blazars of unknown type, 8 misaligned AGNs, 4 RL-NLSyl galaxies, 10 AGNs of other 
types, and 2 starburst galaxies (NGC 4945 has fallen out of the list). 

Important questions that are impacted by or arise from the Fermi-hAT observations of blazars 
and radio galaxies, and interpretation that relates to these questions within the context of 
leptonic models, will occupy the remainder of this short review. 

2. Important questions 

After summarizing the standard leptonic model for blazars, we consider some important 
questions raised by the Fermi-LAT data, including 

(i) Limitations of one- zone synchrotron/SSC model for TeV blazars 

(ii) Relation between aligned and misaligned jet sources 

(iii) Validity and meaning of the blazar sequence 

(iv) GeV spectral cutoffs of LSP and ISP blazars 

(v) Super-Eddington luminosities of blazars 

2.1. Leptonic Models for blazar SEDs 

The general outlines of blazar leptonic models are now well established for the standard one-zone 
model within the context of blazar unification [9]. In this paradigm, the basic ingredients of a 
blazar are an accretion disk, regions of broad and narrow emission-line gas, and a dusty IR torus. 
Radio-loudness is a consequence of collimated relativistic outflows of plasma directed along the 
polar axis of the accretion disk. When observing at small angles to the plane of the accretion 
disk, the torus gas conceals the broad-line region formed close to the black hole, so that only 
narrow emission lines are detected |10j . These are the Type II Seyfert galaxies or, if radio-loud, 
narrow line radio galaxies. Observations at larger angles to the plane of the accretion disk, 
where the torus no longer obscures the central region, reveal the broad emission lines. These 
are the Type I Seyferts or, if radio-loud, the broad-line radio galaxies. 

When viewed at very small angles with respect to the jet axis, the Doppler-boosted 
nonthermal radiation overwhelms the accretion-disk radiation and makes the familiar two- 
humped broadband blazar continuum |11) . The featureless continuum is made by nonthermal 
particles accelerated at shocks formed, for example, by colliding shells [12] or through sheared 
flows in structured jets [13]. In the standard blazar model for flaring activity, a compact plasma 
zone filled with nonthermal electrons and positrons radiate synchrotron photons that form the 
highly polarized radio continuum that extends to optical frequencies in flat spectrum radio 
quasars (FSRQs), and to X-ray frequencies in the high-synchrotron peaked (HSP; with peak 
synchrotron frequency z^f^^ > 10^^ Hz) BL Lac objects. Accompanying the synchrotron emission 
component is the synchrotron self-Compton (SSC) 7-ray component. This simple one-zone 
model successfully fits the SEDs of HSP BL Lac objects like Mrk 421 [14J or Mrk 501 [15j. 
The synchrotron/SSC model for these objects implies Doppler factors (5d ~ 10 - 50, comoving 
magnetic fields B' ~ 0.01 - 0.1 G, and absolute powers ~ 10*^^ erg s^^. The outflow Lorentz 
factors are compatible with values determined by superluminal radio observations (16] . and the 



inferred jet powers are far below the Eddington luminosity for black holes with mass ~ 10 Mq 
that are thought to power these objects. 

The simple one-zone synchrotron/SSC model has difficulty to explain the SEDs of FRSQs 
and intermediate (ISP; 10^"^ Hz < i^p^^ < 10^^ Hz) BL Lacs like BL Lac itself [17j, which 
have dominant 7-ray components that would require, for a synchrotron/SSC model, parameters 
far from equipartition, which would consequently place large demands on jet power. External 
radiation fields, either from the disk [IH], the broad- line region (BLR), or the IR torus [12], are 
then invoked to provide seed photons that are Compton-scattered to form an additional 7-ray 
emission component. Within this framework, precise detailed fits are made to sources like 3C 
279 [TT]. Blazar modeling seemed somewhat satisfactorily explained at the close of the EGRET 
era with synchrotron/SSC models for X-ray selected BL Lac objects, and with the addition 
external Compton scattering components for FSRQs and radio-selected BL Lac objects. 

2.2. Difficulties of simple blazar models 

Even before the launch of Fermi, cracks in our understanding of the blazar phenomenon became 
apparent from observations with TeV telescopes. The July 2006 flaring data of PKS 2155- 
305 observed with HESS [20J and Swift required 60 ^ 100 for good synchrotron/SSC model 
fits [2T]. Very high energy (VHE; > 100 GeV) emission from the FSRQ 3C 279 observed 
with the MAGIC telescope [22] was found to pose severe difficulties to either synchrotron/SSC 
or external Compton scattering models [23]. Unexpectedly low extragalactic background light 
(EBL) intensities at the level of that implied by galaxy counts were required to make sense of the 
unusually hard deabsorbed TeV spectrum of H 2356-309 and lES 1101-232 |24j. Furthermore, 
doubts about the 7-ray emission site being located within the BLR were raised by radio/7-ray 
correlations [25] and multwavelength observations of BL Lacertae, suggesting emission in an 
inner acceleration zone and at an outer recollimation zone J26J. 

These difficulties have been compounded with the advent of Fermi and long baseline 
GeV data correlated with targeted campaigns at radio, optical. X-ray, and VHE energies. 
Broadband modeling [57] of a large sample of TeV BL Lacs found that most are well fit with a 
synchrotron/SSC model with B' ~ 0.01 - 1 G and (5d ~ 20 - 50, but with a second group of 
blazars displaying unusually low magnetic fields and large break electron Lorentz factors. These 
latter sources include lES 1101-232, lES 0229+200, and lES 0347-121, and are peculiar for BL 
Lac objects by being very weakly variable or nonvariable, and displaying very weak GeV fluxes. 
They are now being used to make inferences about the level of the intergalactic magnetic field 
(IGMF) [28j, insofar as the absorbed primary TeV radiation make a large fiux of secondary 
electrons and positrons from attenuated 7 rays. These pairs would make GeV radiation by 
Compton scattering photons of the CMBR and would overproduce the measured Fermi fiux 
unless the pairs are deflected by a sufficiently strong intergalactic magnetic field. The implied 
IGMF ranges from ^ 10" G to 10"^^ G [291 EHl EI] depending on the IGMF correlation 
length and duration of the blazar engine. 

If the IGMF is ~ 10~^^ - 10~^* G, then a separate, quiescent GeV emission component 
will be present in the SED of HSP BL Lac objects. Indeed, this could account for the slowly 
varying GeV fiux in sources like Mrk 501 ^2]. The situation is further complicated if TeV 
blazars accelerate ultra-high energy cosmic rays (UHECRs) and produce a separate high-energy 
radiation component due to photopair production by UHECR protons while traveling through 
intergalactic space [33^ . This could affect inferences of the EBL based on the assumption 
that the highest energy photons are made at the source. This ambiguity can be resolved by 
observations of multi-TeV radiation by CTA and HAWC from sources like lES 0229+200, which 
can be used to discriminate between a photon-induced and UHECR proton-induced origin of 
the high-energy radiation and determine whether TeV blazars are sites of UHECR acceleration 



For FSRQs such as PKS 1510-089 [36] and 3C279 [37], coherent optical polarization position- 
angle swings over timescales of many weeks with late-time flares coinciding with strong 7-ray 
activity seem to be most simply explained if the emission is formed in a helical jet on multi-pc 
scales. It is not clear, however, how this relates to indications in PKS 1510-089 that the 7-ray 
flux precedes the optical flux by ~ 13 d [38j. Moreover, the interpretation of distant emission 
regions may collide with the inference about sharp GeV cutoffs in the SEDs of these sources, as 
described in the next section. 

2. 3. Ge V spectral cutoffs of LSP and ISP blazars 

Already by the time of the LBAS and the first paper on 3C 454.3 [39] it was evident that the 7- 
ray spectra of FSRQs, which are nearly all low synchrotron-peaked (LSP) sources with uFi, peak 
synchrotron frequency u^^^ < 10^^ Hz, as well as ISP objects such as AO 0235-1-164, display 
remarkably sharp cutoffs between 1 and 10 GeV. This feature is particularly evident in the 
compilation of spectra of Fermi bright blazars @0j, where a distribution of source- frame break 
energies as a function of apparent isotropic 7-ray luminosity is displayed. Essentially 
all sources have E^^ within 3a of 5 GeV. Two other important aspects of this result should be 
mentioned: (1) The break energy changes only slightly with flux state. In the case of 3C 454.3, 
E^^ increases by a factor of 2 when the flux varies by more than a factor of 20, as indicated by 
spectral analysis of the bright GeV blazar 3C 454.3 [5l]. (2) The change in spectral index below 
and above the break is typically larger than 1 unit, ruling out naive radiative cooling models 
that predict a spectral break of 0.5 units. 

This has generated considerable theoretical interest. A hybrid scattering scenario |42j 
involving both accretion-disk and BLR photons can reproduce the SED of 3C454.3, with the 
added feature that the spectral break energy is insensitive to location if the BLR has a wind-like 
profile, because the energy densities of both the accretion-disk and BLR radiation fields then 
scale with the same dependence on radius. Poutanen and Stern argue in print [43 ^ .44j and at 
this conference that He II recombination line and continuum radiation attenuate a continuum 
GeV emission component through 77 pair production. For He II recombination radiation at 
energy sa 4 x 13.6 = 54.4 eV or = 54.4/511000 = 1.06 x 10""^ (in mgC^ units), photons with 
energy mgC^ei > 2/e = 9.6 GeV, close to E^^, can be strongly absorbed. The 7-ray emitting 
region must be found deep within the BLR for this model to work, and the excessive target 
photon energy density may introduce difficulties to fit the multi-band spectral data of 3C 454.3. 

A further possibility is that the breaks are due to Klein-Nishina effects on Compton-scattered 
Ly a photons [H]. GALEX observations of 3C 454.3, as reported by Bonnoli and coworkers 
|45j . reveal the presence of a strong Ly a radiation field in the BLR of 3C 454.3, amounting 
to a line luminosity of (2 - 4)xl0^^ erg s~^. The onset of the Klein-Nishina decline takes 
place when A'je = 1 for electrons with Lorentz factor 7 scattering Ly a target photons with 
dimensionless energy e = 10.2/511000 = 2 x 10^^. Photons are Compton-scattered to energies 

= (3/2)7^ e, implying ee^ = 3/32, i.e., m^c^es = 2.4 GeV, consistent with the measured break 
energy. Spectral calculations employing a power-law electron spectrum produce a break at the 
right energy, but with an insufficiently sharp break [51] • Combined Compton and synchrotron 
losses will harden the electron spectrum and reduce the magnitude of the break [4:6\ I47j . but a 
curved electron spectrum described by a log parabola function [18], formed when acceleration 
dominates over cooling, might resolve this difficulty. 

2.4- Blazar sequence 

The search for an organizing principle in blazar studies has centered on the so-called blazar 
sequence, which is the well-known anti-correlation of the I'Fi, synchrotron and Compton peak 
frequencies with luminosity [50\ I51j . It is reflected in the Fermi data by plotting the blazar 
spectral index F^ versus the 7-ray luminosity. This figure reveals a clear separation between 



Figure 1. Gamma-ray spec- 
tral slopes of BL Lac 
objects (open blue circles), 
FSRQs (open black squares), 
FRl radio galaxies (red cir- 
cles), FR2 radio sources (green 
squares) , and star- forming 
galaxies (magenta diamonds), 
are plotted as a function of 
their 100 MeV - 5 GeV 7- 
ray luminosity [39]. The 7- 
ray emitting misaligned AGNs 
are ILAC sources associated 
with objects in the Third 
Cambridge Revised (3CRR) 
and Molonglo radio catalogs. 
The 7-ray emission from star- 
forming galaxies are powered 
by cosmic-ray processes rather 
than AGNs. 

BL Lac objects and FSRQs — a blazar divide |52J — that is attributed to a difference in accretion 
regime. As plotted in Figured] one sees that all FSRQs have F^ > 2. Inasmuch as the presence 
of strong emission lines means that there has to be a significant surrounding scattered radiation 
field, then a consistent interpretation of FSRQ SEDs is to attribute the 7-ray emission to external 
Compton processes [5T1, with the trend in the SEDs from FSRQ to BL Lac objects resulting from 
the reduction in the intensity of the surrounding radiation field [5l3J • This is reasonable provided 
that the emission region of FSRQs is found within the BLR. Indeed, this seems to be the most 
natural explanation for rapid variability [53]. Yet this expectation collides with detection of 
VHE 7 rays from FSRQs like 3C 279 ^22j and, most recently and most spectacularly, 4C +21.35, 
which exhibits a 10 minute flaring episode of 70 - 400 GeV radiation [SS] from a source with 
very strong IR radiation |56j that would strongly attenuate high-energy radiation made within 
the pc scale. Two-zone models have been advocated to deal with these questions |23t 157) . 

The blazar sequence and blazar divide are not, however, immune to criticisms, most notably 
that a large fraction of BL Lac objects lack redshifts. This fraction amounts to more than 50% 
in the 2LAC. Because strongly beamed nonthermal radiation would conceal emission lines, then 
the highest redshift, highest luminosity BL Lac objects would be selectively excluded from the 
blazar sequence diagram, and these would also likely be the HSP BL Lacs where the nonthermal 
radiation is strongest in the optical/UV range. Besides outliers to the blazar sequence, there is 
evidence for a population of HSP FSRQs [58j. A further concern is that the BL Lac objects are 
primarily identified from and biased towards the radio and X-ray surveys from which they are 
selected [59] , 

In any case, the simplest characterization of the blazar sequence must consider the 
contribution of slightly off-axis jetted sources and misaligned radio galaxies. When observing 
away from the beaming axis, the apparent luminosities from the smaller Doppler boosting should 
appear at lower peak frequencies, opposite to the blazar sequence correlation. Consequently 
radio galaxies like NGC 6251 [60] will fill in a distinct portion of a diagram of the blazar 
sequence, turning it rather into an envelope [61j . 
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2. 5. Relation between aligned and misaligned jet sources 

A curiosity revealed by Fig. [1] has to do with the spread in luminosities of the radio galaxies 
compared with the aligned blazar jet sources. The misaligned AGN (MAGN) population 
detected with Fermi consists of 11 sources |49J. Theses includes 7 FRI and 2 FRII radio galaxies, 
which are thought to be the parent population of BL Lac objects and FSRQs, respectively, in 
addition to 2 steep spectrum radio quasars that are believed to be slightly misaligned FSRQs. 
The MAGN sources in the 3CRR catalog have large core dominance parameters compared to the 
general 3CRR source population, implying that the beamed component makes an appreciable 
contribution to the 7-ray flux. This is furthermore supported by the fact that three of these 
sources — 3C 78, 3C 111, and 3C 120 — do not appear in 2LAC, evidently due to variability. 

Why the ratio of measured 7-ray luminosities of FRI galaxies and BL Lac objects span a 
much larger range than the comparable ratio for FRII galaxies and FSRQs is an open question, 
if in fact it is not due to statistics. The simplest possibility is that there are a lot more nearby 
FRIs, and limitations on Fermi sensitivity will therefore favor detection of these nearby sources. 
At 2 -C 1, the luminosity distance di ~ 4200^(1 + z) Mpc. The Fermi-LAT reaches a limiting 
energy-flux sensitivity of ~ 5 x 10^^^ erg cm^^ s^^ for two years of observations that is, unlike 
integral photon flux, only weakly dependent on source spectral index [8]. Fermi can thus only 
detect sources with luminosity ^ 10^^2:^(1 + 2z) erg s~^. The number of nearby FRI and 
FRII galaxies depends on the space density of these objects which, in the study by Gendre, 
Best, and Wall [62], is inferred from 1.4 GHz NVSS-FIRST radio observations (recall that the 
Fanaroff-Riley dichotomy is based on 178 MHz luminosities P178 mHz)- From their Figure 12, 
FRI and FRII radio galaxies have a local {z < 0.3) space density of w 35 x 10~^ Mpc~^ and 

2 X 10^^ Mpc~^, respectively, at Pi. 4 gHz ~ 10^^ W Hz~^ sr^^, implying that the local space 
density of FRIs exceeds that of FRIIs by ~ 20. 

This implies that there are several hundred FRIs and FRIIs within z = 0.1, yet other than 
3C 111, all LAT-detected radio galaxies within this volume are FRI sources (see Figure 3 of 
[63]). This is further remarkable in that FSRQs and their putative parent population of FRII 
radio galaxies are far more radio luminous than BL Lac objects and FRI galaxies. With only 
a few hundred randomly aligned sources within z = 0.1, it might be the case that a narrower 
7-ray beaming cone in FSRQs, with a more rapid fall-off in off-axis flux, makes detection of 
these nearby sources far less likely than a broader 7-ray emission cone in BL Lac objects, as 
expected for the different beaming factors from SSC emission and external Compton processes 
[Ml EH]. This could also reflect differences in jet structure between FSRQs and BL Lac objects 
[66], or extended jet or lobe emission in FRIs that is missing in FRII galaxies. 

2.6. Apparent super- Eddington luminosities of blazars 

We end this contribution on a speculative note. The time-averaged luminosities of FSRQs extend 
to values in excess of ~ 10^^ erg s^^ (Figure [1]). In the extraordinary 2010 November flare of 
3C 454.3 [67], reached apparent isotropic luminosity of (2.1 ±0.2) x 10^*^ erg over a period 
of a few hours, making it the most luminous blazar yet observed. Black-hole mass estimates 
for 3C 454.3 are in the range 0.5 < Mg < 4 [ig, where lO^MgM© is the mass of the black hole 
powering this AGN. For this range of masses, the Eddington luminosity therefore ranges from 
~ 6 X 10^^ erg s~^ to ~ 5 x 10^^ erg s~^. During this extreme outburst, the apparent luminosity 
of 3C 454.3 was more than a factor of ~ 400 greater than its Eddington luminosity. Even its 
time-averaged luminosity of ~ 10^^ erg s~^ is super-Eddington by a factor of ~ 20. 

Assuming that the Eddington condition does limit accretion flow onto the black hole, which 
is likely to be the case for the long-term average luminosity if not for the flaring luminosity, then 
the absolute radiant luminosity is limited to a value of Labs ~ 5 x 10^^ erg s"^. This is consistent 
with the large apparent luminosities if the emission is highly beamed. For a simple top-hat jet 
beaming factor, a jet opening angle 6j implies a beaming factor fh = G'j/2 for a two-sided jet 



with 6j ^ 1. A mechanism for colhmation is, however, required. Should this arise from the 
Blandford-Znajek process, then we are still restricted to values of the absolute Blandford-Znajek 
power 

^lBZ-vrc(-) r^i?^ ' 

where a/M is the angular momentum per unit mass, is the radius of the event horizon, and 
-B!j_ is the radial component of the magnetic field threading the event horizon [BB]- Scaling the 
energy density of the magnetic field to the energy density of accreted matter near the event 
horizon shows that the Blandford-Znajek power is likewise Eddington-limited. Making the 
hypothesis that the extraction of energy through black-hole rotation collimates the jet outflow 
with cos6'j ^ a/M, then fb ^ l-{a/M) and 9j ^ y/2{l - a/M), implying a/M > l-{LEdd/L^)- 
If the jet opening angle is a consequence of the bulk Lorentz factor T of the outflow, then 

r ^ ^ L^/2LEdd- 

For the case of 3C 454.3 in its flaring state, when L^/LEdd ~ 10^, this hypothesis then implies 
that a/M > 0.999 and F ^ 23, which is consistent with the value F^m ~ 14 from 77 opacity 
arguments |67) . In the most conservative case with Mg = 4, a/M = 0.998 and F ^ 15. This is 
marginally consistent with the limiting maximum value a/M = 0.996 suggested by Aschenbach 
from analyses of microquasars and the Galactic Center black hole [69]. The smaller black- hole 
mass estimate, Mg = 0.5 (Bonnoli et al. 2010), implies a value of a/M that violates this limit by 
a large margin. Much work, both numerical and theoretical, has been devoted to jet formation 
from the Blandford-Znajek process, and it is unclear if jet collimation can be described by the 
guess that cosOj a/M, but it is interesting to suggest this possibility, which leads to values of 
F consistent with separate inferences regarding the outflow Lorentz factor. 
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